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ABSTRACT 

General relativity predicts that the spin axes of the pulsars in the double- pulsar system 
(PSR J0737— 3039A/B) will precess rapidly, in general leading to a change in the observed pulse 
profiles. We have observed this system over a one-year interval using the Parkes 64-m radio 
telescope at three frequencies: 680, 1390 and 3030 MHz. These data, combined with the short 
survey observation made two years earlier, show no evidence for significant changes in the pulse 
profile of PSR J0737— 3039A, the 22-ms pulsar. The limit on variations of the profile 10% width 
is about 0?5 per year. These results imply an angle 5 between the pulsar spin axis and the orbit 
normal of < 60°, consistent with recent evolutionary studies of the system. Although a wide 
range of system parameters remain consistent with the data, the model proposed by Jenet & 
Ransom (2004) can be ruled out. A non-zero ellipticity for the radiation beam gives slightly but 
not significantly improved fits to the data, so that a circular beam describes the data equally 
well within the uncertainties. 

Subject headings: pulsars: general — pulsars: individual (J0737-3039A) 



1. Introduction 

The double-pulsar system PSR J0737-3039A/B 
(Burgay et al. 2003; Lyne et al. 2004) provides a 
wonderful laboratory for investigations of rela- 
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tivistic gravity (Kramer et al. 2004) as well as the 
physics of pulsar magnetospheres (e.g., McLaugh- 
lin et al. 2004). The system consists of a 22-ms 
pulsar, PSR J0737-3039A, in a 2.4-h binary orbit 
with PSR J0737— 3039B, a younger pulsar with a 
spin period of 2.7 s. The system is mildly eccentric 
(e ~ 0.088) and is viewed nearly edge-on (orbit 
inclination i ~ 88°). With mean orbital speeds 
v ~ 0.001c, the system is highly relativistic, al- 
lowing the detection of four "post-Keplerian" pa- 
rameters in just six months of observation (Lyne 
et al. 2004). The post-Keplerian parameters, to- 
gether with the mass ratio, uniquely measurable 
in this double-pulsar system, give accurate values 
for the masses of the two stars as well as stringent 
tests of general relativity (Kramer et al. 2004). 

Within the framework of general relativity, the 
spin vectors of the two pulsars are expected to 
exhibit geodetic precession about the total angu- 
lar momentum of the system (Damour & Ruffini 
1974; Barker & O'Connell 1975). Since the to- 
tal angular momentum is dominated by the or- 



bital motion, the pulsar spin vectors effectively 
precess about the orbit normal. The precession 
rate depends on the Keplerian parameters and 
the masses of both pulsars, with predicted pre- 
cessional periods of ~ 75 yr and ~ 71 yr for 
PSR J0737-3039A and PSR J0737-3039B re- 
spectively (Lyne et al. 2004). These are about 
a factor of four shorter than the value for PSR 
B1913+16, the Hulse- Taylor binary pulsar (Wcis- 
berg, Romani & Taylor 1989). In general, pre- 
cession will result in a change of the width and 
shape of the observed profile as the angle of the 
line of sight with respect to the pulsar spin axis 
(£) changes during the precessional period. Long- 
term evolution of the mean pulse profile has been 
observed in PSR B1913+16 (Weisberg et al. 1989; 
Kramer 1998) and PSR B1534+12 (Stairs et al. 
2004) and interpreted as evidence for geodetic pre- 
cession. This gives the three-dimensional geome- 
try (within a 180° ambiguity) of the system and, 
in the case of PSR B1913+16, a prediction that 
the pulsar will no longer be visible after the year 
2025. Weisberg & Taylor (2002) have used the 
precessional motion to produce a two-dimensional 
image of the part of the PSR B1913+16 emission 
beam traversed so far, suggesting that the beam 
is elongated in the latitudinal direction. 

If the magnetic-pole model (Radhakrishnan 
& Cooke 1969) is assumed, ( and the inclina- 
tion angle of the magnetic axis with respect to 
the spin axis (a) can be deduced from a fit of 
the model to the observed variations of polar- 
ization position-angle, at least for wide profiles 
(e.g., Lyne & Manchester 1988). Demorest et al. 
(2004) used such an analysis to conclude that 
PSR J0737-3039A is a nearly aligned system with 
a <~ 4°, but the fit to the observed position angle 
variations is poor. The observer angle £ is uncon- 
strained by this fit. Jenet & Ransom (2004) deter- 
mined possible geometries of the J0737— 3039A/B 
system based on a model for the orbital mod- 
ulation of the PSR J0737— 3039B pulse inten- 
sity. Two solutions were obtained with (a, 5) of 
(1?6 ± 1?3, 167° ± 10°) and (14° ± 2°, 90° ± 10°) 
respectively. The first of these is consistent with 
the Demorest et al. (2004) solution. Both so- 
lutions predict a rapid evolution of the observed 
profile width of PSR J0737-3039A, with expected 
changes of ~ 42° and ~ 96° per year for the two 
models. Here we report on observations of the 



mean pulse profile of PSR J0737— 3039A over a 
three-year interval, 2001 August to 2004 August. 

2. Observations and Analysis 

We have used the Parkes 64-m radio telescope 
to observe the PSR J0737-3039A/B system since 
2003 May 1 (MJD 52670), shortly after its confir- 
mation, to 2004 August 8 (MJD 53225). Observa- 
tions were made in three frequency bands centered 
at 680, 1390 and 3030 MHz, with the 680-MHz 
(50cm) and 3030-MHz (10cm) observations com- 
mencing in 2003, December. Before 2003 October, 
the 1390 MHz observations were made using the 
center beam of the multibeam receiver, which has 
an equivalent system noise of approximately 29 Jy. 
In addition, we have the original survey pointing, 
a 4-min observation made at 1390 MHz on 2001 
August 22 (Burgay et al. 2003). After 2003 Oc- 
tober, the H-OH receiver with a system noise of 
about 42 Jy was used. Both the multibeam and 
H-OH systems were dual polarization and used a 
filterbank system consisting of2x512x0.5 MHz 
channels. After detection, the two polarizations 
for each channel were summed, sampled at either 
80 lis or 125 /is intervals and recorded to tape 
for subsequent processing. The other two bands 
were observed using the dual-polarization coaxial 
10cm/50cm receiver, with system noises of <~ 48 
Jy and ~ 64 Jy at 10cm and 50cm respectively. A 
filterbank system with 2 x 256 x 0.25 MHz channels 
was used with the 50cm receiver and at 10cm the 
filterbank system had 2 x 192 x 3 MHz channels. 

Over the 15-month interval, we obtained 59 ob- 
servations at 1390 MHz, 17 observations at 680 
MHz and 27 at 3030 MHz, all with durations 
ranging from 10 min to 5 hr. Data from each 
observation were folded at the apparent topocen- 
tric period of PSR J0737-3039A to form mean 
pulse profiles with 256 phase bins, which were then 
summed in frequency and time to form a single 
profile using the PSRCHIVE data analysis sys- 
tem (Hotan, van Straten & Manchester 2004). All 
observations for each frequency band were then 
summed to form grand-average profiles containing 
116, 27 and 34 hrs of data, respectively. These 
profiles are shown in Figure 1. 

The 1390-MHz observations over the 15 months 
from 2003 May were divided into eight chronologically- 
ordered groups and average profiles were formed 
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Fig. 1. — Mean pulse profiles for 

PSR J0737-3039A at three frequencies, where 
360° of longitude corresponds to the pulse period. 
The 680-MHz and 3030-MHz profiles were aligned 
to give maximum cross-correlation with the 
1390-MHz profile. 

for each of these eight groups. Integration times 
ranged between 5 hrs (MJD 53068) and 20 hrs 
(MJD 53000). Examination of these average pro- 
files showed little or no evidence for evolution of 
the mean pulse profile over the 15-month inter- 
val. The eight profiles were each aligned with 
the grand-average profile using a cross-correlation 
analysis and baselines subtracted, taking the lon- 
gitude range 310° to 50° (Figure 1) as the baseline 
region. Each profile was scaled to have the same 
area as the grand-average profile, which was then 
subtracted from it to form the difference profiles 
shown in Figure 2. Also shown at the bottom of 
the figure is the difference profile corresponding 
to a stretched version of the grand-average profile, 
demonstrating the signature of a change in profile 
width. The adopted stretch factor (1.004) corre- 
sponds to an increase of about 1° in profile width, 
where the width is defined to be the separation 
of the steep outer edges of the profile (at roughly 
60° and 300° in Figures 1 and 2) at 10% of the 
amplitude of the second (stronger) component. 

Significant differences are seen, especially be- 
tween MJDs of 53000 and 53103. In general, these 
differences do not have the signature of an increase 
or decrease in profile width, but rather appear to 
result from changes in the relative amplitude of the 



various pulse components, especially at the lead- 
ing and trailing components of the profile. Polar- 
ization observations (e.g., Demorest et al. 2004) 
show that these components have high linear po- 
larization. It is likely that the observed changes 
are an instrumental effect resulting from differ- 
ences in receiver gain between the two polariza- 
tion channels combined with parallactic angle vari- 
ations. 

These difference profiles confirm the lack of sig- 
nificant secular profile evolution. There is no evi- 
dence for a systematic change in the overall profile 
width to a level of much less than 1° and there is no 
evidence for any intrinsic changes in the width or 
shape of the leading and trailing pulse components 
considered separately. We have also reprocessed 
the original 4-min survey pointing to give a pro- 
file with signal-to-noise ratio of about 25. Again, 
this shows no evidence for any significant change 
in profile shape or width. 

In order to quantify the profile stability, widths 
at 10% of the pulse peak were measured for each 
of the eight 1390-MHz profiles and the 2001 ob- 
servation using a cross-correlation technique. The 
leading (50° - 130° in Figures 1 and 2) and trail- 
ing (240° - 310°) pulse components were sep- 
arately cross-correlated with a standard profile 
derived from the grand-average profile contain- 
ing just the relevant component, determining the 
phase of each component relative to a reference 
phase. The difference of these phases was then 
added to the 10% width of the grand-average pro- 
file to give the 10% width at each epoch. Errors 
from the cross-correlation analysis were increased 
in quadrature by 10 /is to allow for systematic ef- 
fects. The derived 10% widths are plotted in Fig- 
ure 3 along with the fitted trend line. Clearly, 
there is no significant trend, with the fitted slope 
being -0?2 ± 0?3 per year. The 10% profile width 
at MJD 53000 (2003 December 27) from the fit is 
238?45±0?12. Quoted errors arc ±2<r for both the 
width and slope. 

3. Discussion and Conclusions 

Figure 1 confirms the overall bilateral symme- 
try of the PSR J0737-3039A pulse profile, consis- 
tent with its origin from magnetic field lines asso- 
ciated with a single magnetic pole (Burgay et al. 
2003; Demorest et al. 2004). There are significant 
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Fig. 2. — Difference profiles for 

PSR J0737-3039A at 1390 MHz relative to 
the grand-average pulse profile shown at the top 
of the figure for eight epochs, together with the 
difference profile corresponding to the grand- 
average profile stretched by 0.4%. The vertical 
scale of the difference profiles is a factor of ten 
larger than that of the grand-average profile. 

differences in the spectral indices of the various 
pulse components. The most striking of these is 
the relatively flat spectrum of the component at 
the extreme leading edge of the profile which is 
relatively much stronger at 3030 MHz. The 3030- 
MHz profile also shows a matching flat-spectrum 
component at the extreme trailing edge of the pro- 
file, reinforcing the idea that both components of 
the observed profile originate on field lines associ- 
ated with a single magnetic pole. We assume this 
in the following discussion. 

The main result of this paper is the extraor- 
dinary stability of the PSR J0737-3039A pulse 
profile over more than three years, a significant 
fraction of the expected geodetic precessional pe- 
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Fig. 3. — Profile widths at 10% of the peak am- 
plitude for PSR J0737-3039A at 1390 MHz as a 
function of MJD relative to 53000. Error bars are 
±lcr. The dashed line is a formal weighted least- 
squares fit of a linear trend to these data. 



riod of 75 years. Provided the angle between the 
spin axis of the pulsar and the orbital angular mo- 
mentum (5) is non-zero, we would expect to see a 
variation in profile width as a function of time. 
In particular, the system geometries proposed by 
Jenct & Ransom (2004) predict changes of many 
tens of degrees per year in the observed profile 
width. Clearly, these are not observed, ruling out 
this model in its present form. 

One possible explanation for the observed lack 
of variation in profile width is that the spin axis 
of PSR J0737-3039A is aligned with the orbital 
angular momentum, i.e., 5^0. Since we view 
the orbit nearly edge-on, even for an orthogonal 
magnetic axis (a ~ 90°), this would require an ef- 
fective beam radius p > 90°, which seems unlikely. 
Furthermore, given that PSR J0737-3039B prob- 
ably suffered a significant natal kick, an aligned ro- 
tation axis for PSR J0737-3039A is possible but 
unlikely (Willems, Kalogera & Henninger 2004). 

We have investigated the limits which can be 
placed on the geometry of the system by modeling 
the observed profile width as a function of time 
(cf. Kramer 1998). An inclination angle i = 88° is 
assumed 1 . For a circular emission beam, four pa- 
rameters are sufficient to model the expected pro- 
file changes: a, 8, p and an epoch T describing the 



1 Since only sini is determined at present, equivalent solu- 
tions arc possible with i = 92° and a — > 180° — a, etc. 
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precession phase. We fitted a circular beam model 
with p < 90° to the data and derive a x 2 -sphere 
plot for a and 5 as shown in Figure 4. This fig- 
ure shows that solutions with misalignment angle 
S < 60° are preferred, consistent with the conclu- 
sions of Willems et al. (2004) based on the ob- 
served velocity of the system. 

The formal best-fit solution (reduced x 2 = 
2.88) is located at a ~ 19° and 6 ~ 14° (and 
at a corresponding mirror-solution of a ~ 161° 
and 8 ~ 166°). The predicted time variation for 
the observed profile width for this configuration 
is shown in Figure 4(a); pulses would be detected 
over the whole precessional period with an approx- 
imately sinusodial variation in the profile width. 
However, the x 2 -sphere is rather flat, and a statis- 
tically satisfactory solution is possible over a wide 
range of angles. Figure 4(b) shows a solution with 
small a and S which is consistent with the pre- 
ferred solution of Dcmorest et al. (2004) based on 
the observed position angle variations. Both of 
these solutions have a beam radius of very close 
to 90° corresponding to a fan beam. In case (b) 
however, the pulse is visible for a limited period 
only roughly 1984 to 2023. Figure 4(c) shows a 
solution with intermediate a and S for which the 
beam radius is smaller, about 65°. In all solutions, 
the observed pulse is emitted from field lines asso- 
ciated with a single magnetic pole as indicated by 
the observed profile symmetry 

Unless S ~ 0, the essentially constant profile 
width observed so far implies that the system is 
at or close to precessional phase zero (or 180°). 
At these times is dW/dt ~ 0. 2 This is possi- 
ble - a similar situation evidently applies for PSR 
B1913+16 (Kramer 1998) - but statistically un- 
likely. A possible way of relaxing this constraint 
is to allow a non-circular beam. Biggs (1990) and 
Kapoor & Shukre (1998) have shown that, under 
certain assumptions, the polar emission beam is 
compressed in the latitudinal direction for a > 0°. 
On the other hand, polarization observations of 
young and millisecond pulsars (e.g., Narayan & 
Vivekanand 1983; Manchester & Han 2004) sug- 
gest beams effectively elongated in the latitude 
direction. First, we took the model of Kapoor & 

2 There is a special case with cos p sin S cos <t> = cos a, where 
<E> is the precessional phase, which results in dW/dt = 0, 
but this requires an improbable fine-tuning of the system 
parameters. 



Shukre (1998) in which the beam compression is a 
function of a only and recomputed the x 2 -sphere. 
The results were only marginally different to those 
shown in Figure 4. We then took the beam ellip- 
ticity as a free parameter, first doing a grid search 
in ellipticity at each point on the (a, 5) plane and 
then searching for a global x 2 minimum in (p m , 
To, e) space, where p m is the beam major axis, To 
is the time of precessional phase zero and e is the 
beam ellipticity. Latitudinally compressed beams 
with axial ratio piat/piong ^ 0.5 are restricted to a 
fairly small range of a and 5 around 45°. On the 
other hand, fan-like beams with large axial ratio 
allow a much larger range of a and S with values 
< 30° preferred. However, in none of these cases 
was the x 2 value significantly better than those 
found in the circular-beam case. 

In conclusion, we find that the pulse-width 
variations (or lack of them) observed so far al- 
low a wide range of configurations for the PSR 
J0737-3039A/B system. Models with non- 
circular beams give somewhat better fits to the 
data but are statistically indistinguishable from 
fits with a circular beam. Misalignment angles 
S < 60° are generally preferred, consistent with 
the conclusions of Willems et al. (2004), and 
the configurations discussed by Jenet & Ransom 
(2004) can be ruled out. Solutions in which the 
pulse never disappears from view are possible, 
as are solutions where it does, in the most ex- 
treme cases about 15 years from now. Clearly, 
a longer time baseline, even one more year, will 
help to constrain the models, as will a more real- 
istic interpretation of the observed position-angle 
variations. 
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Fig. 4. — Results of fits to the observed profile 
widths of a model for geodetic precession of the 
pulsar spin axis, inclined by an angle 5 to the or- 
bit normal, with a circular beam of radius p in- 
clined at an angle a to the pulsar spin axis. The 
upper part shows contours of \ 2 a t 1-f intervals, 
relative to the minimum value, with the darker re- 
gions having lower \ 2 ■ The lower cutoff to the al- 
lowed regions represents the locus of points where 
the beam radius p = 90°. These contours are mir- 
rored for 90° < 5 < 180° with (a, S) corresponding 
to (180° - a, 180° - 5). In the lower part of the 
figure, the variation of profile width is shown as a 
function of time for three cases: (a) a point near 
the global minimum for x 2 , (b) a small a, small 5 
case, and (c) a solution with intermediate a and 
5. The angle ( is the minimum value of the incli- 
nation of the line of sight to the pulsar spin axis 
which occurs at precessional phase zero (or 180°) 
when the pulsar spin axis is in the plane defined 
by the orbit normal and the line of sight. 
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